The Poisson–Boltzmann equation is a differential equation that describes electrostatic interactions between molecules in ionic solutions. It is the mathematical base for the Gouy–Chapman double layer (interfacial) theory; first proposed by Gouy in 1910 and complemented by Chapman in 1913. The equation is important in the fields of molecular dynamics and biophysics because it can be used in modeling implicit solvation, an approximation of the effects of solvent on the structures and interactions of proteins, DNA, RNA, and other molecules in solutions of different ionic strength. It is often difficult to solve the Poisson–Boltzmann equation for complex systems, but several computer programs have been created to solve it numerically.
The equation can be written as (in cgs):
where is the divergence operator, represents the position-dependent dielectric, represents the gradient of the electrostatic potential, represents the charge density of the solute, represents the concentration of the ion i at a distance of infinity from the solute, is the charge of the ion, q is the charge of a proton, is the Boltzmann constant, T is the temperature, and is a factor for the position-dependent accessibility of position r to the ions in solution. If the potential is not large compared to kT, the equation can be linearized to be solved more efficiently, leading to the Debye–Hückel equation.[1][2][3]